Geometrical properties of Potts model during the coarsening regime 
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We study the dynamic evolution of geometric structures in a poly-degenerate system represented 
by a ^-state Potts model with non-conserved order parameter that is quenched from its disordered 
into its ordered phase. The numerical results obtained with Monte Carlo simulations show a strong 
relation between the statistical properties of hull perimeters in the initial state and during coarsening: 
the statistics and morphology of the structures that are larger than the averaged ones are those of 
the initial state while the ones of small structures are determined by the curvature driven dynamic 
process. We link the hull properties to the ones of the areas they enclose. We analyze the linear 
von-Neumann-Mullins law, both for individual domains and on the average, concluding that its 
validity, for the later case, is limited to domains with number of sides around 6, while presenting 
stronger violations in the former case. 
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I. INTRODUCTION 

After a rapid quench across a phase transition, two or 
more equilibrium states may compete for the growth of 
local order; this reflects in the out of equilibrium evolu- 
tion observed in many different macroscopic systems once 
they reach the dynamic scaling regime [l|-Q. In these 
systems, the properties of macroscopic observables, such 
as correlation functions and linear responses, can still be 
described in relatively simple terms by resorting to the 
dynamic scaling hypothesis. This hypothesis states that 
averaged macroscopic observables depend upon time only 
through a monotonically growing, time-dependent char- 
acteristic length, R{t). This length, whose growth law 
depends on the case under consideration, characterizes 
the average linear size of equilibrated patches at each 
instant t. For example, after being cooled below their 
critical temperature magnetic systems form domains in 
which spins are strongly correlated and the magnetiza- 
tion is uniform. More precisely, on a lattice, a (geomet- 
ric) domain is defined as the ensemble of nearest neighbor 
sites whose spins are aligned. Each of these domains has 
a hull or external border, whose length is the number of 
external sites that are first neighbors of the chosen do- 
main [El and that, in the continuum limit, is referred as 
the perimeter. Figure [1] shows a sketch that illustrates, 
on a square lattice, this and other related definitions. 

Systems as diverse as soap froths [6|-[8|, cellular tis- 
sues and other natural tilings p, [Tq| , magnetic grains 
or poly crystalline structures |lll-[l5| , type-I superconduc- 
tors [16|, etc. are, in a statistical sense, geometrically 
similar. Their overall morphology and growth properties 
are well described by simple spin models with multiple 
ground states. An example is the ferromagnetic Potts 
model [13 whose variables are spins that take q possi- 
ble values and interact with their nearest neighbors. The 
coupling exchanges are taken to be isotropic and homo- 
geneous, i.e. independent of the spin variables and lat- 
tice orientation, favoring spin alignment at low enough 



temperature and, therefore, the existence of q degenerate 
equilibrium states below the phase transition. The criti- 
cal temperature of the bidimensional model on the square 
lattice is known for all values of kj^Tc = 2 J/ ln(l + y^), 
but the nature of the transition depends on being con- 
tinuous and second-order for g < 4 while discontinuous 
for > 4. When submitted to a sudden quench from 
the high to the low temperature phase, and let subse- 
quently evolve with stochastic non-conserved order pa- 
rameter (the magnetization) dynamics, the system tends 
to organize in progressively larger ordered structures of 
each equilibrium state |18l-[22]. One of the most com- 
mon applications of this model is to crystal grain growth, 
where each spin configuration represents a different grain 
orientation and for g ^ oo the continuum limit is real- 
ized 



The simple presence of domains implies the existence 
of topological defects, in this case, the interfaces. When 
the final quench temperature is not too close to zero (to 
avoid pinning effects when q > 2) and much lower than 
the critical temperature (so that the starting ordering 
process does not occur by nucleation when the transi- 
tion is of first order) [25|, [26| , at sufficiently long times 
the walls around large domains tend to be flat and the 
evolution (with non-conserved order parameter) is driven 
by their curvature. Indeed, the surface tension implies a 
force per unit area that is proportional to the local mean 
curvature, acting on each point on the wall. This force 
is responsible for the motion of the interfaces. The field 
theoretical description (a la time-dependent Ginzburg- 
Landau equation) in the continuum leads to the Allen- 
Cahn equation for the local velocity, v = —{Xq/27r)K.ns^ 
where Xq is a temperature and g-dependent dimensional 
constant related to the surface tension and mobility of a 
domain wall and fis is a unit vector normal to the sur- 
face [13 . The sign indicates that the velocity points in 
the direction that tends to reduce the local curvature. 
Thermal effects play a minor role, affecting only the con- 
stant Xq. The time-dependence of the area contained 
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FIG. 1. (Color online.) Illustration of various interfacial 
length definitions on the square lattice. The colored sites 
are all occupied by the same spin species while spins on white 
squares are in a different state. The hull length (solid cir- 
cles) is the external border of the domain. The hull perimeter 
(open circles) contains only the external nearest neighbors of 
the domain (not the internal ones). The domain perimeter 
(crosses) is the ensemble of nearest neighbors of the domain 
and includes the internal ones as well. The measurements 
presented here are for the hull perimeter only. 



within any finite interface can be deduced by integrating 
the velocity around the hull, i.e. the external interface. 
This calculation is specially simple in d = 2 since one can 
use the Gauss-Bonnet theorem to find a close expression 
for the area rate of change, that can be written in the 
following compact form: 



dA 
~dt 



y /\dl = — (n- 
5 



6) if g > 2 , 



(1) 



where n is the number of sides of the external interface, 
each side being the border between the chosen area A 
and one of its neighboring domains. We redefined for 
convenience. While all non-percolating hull-enclosed ar- 
eas of a system with q = 2 states have only one side, 
for g > 2 hull-enclosed areas with more than one side 
can exist and n can be larger than one. In the latter 
case, the equation yields the von Neumann-Mullins law 
for the area A of an n-sided hull-enclosed area [28, 29]. 
As a result, this area can grow, shrink or remain constant 
depending upon the number of sides being larger than, 
smaller than or equal to 6, respectively. Moreover, in the 
course of the evolution, the number of sides that a given 
external interface has can change, so that the equation 
ruling the area evolution changes through the time de- 
pendence in n, a function that one cannot characterize 
in full detail. It is also clear that for g > 2 areas do not 
evolve independently of each other, as occurs for g = 2, 
as n is a quantity that also affects the behavior of the 
neighboring areas. 

For (7 = 2, i.e. the bi-dimensional Ising model, the time- 
dependence of a given hull-enclosed area can be easily 
determined from the integration of Eq. ([1]) since n{t) = 1 
for all areas and these evolve independently. The num- 
ber of hulls with enclosed area A per unit system area 



{A^ Aq, where Aq is a small area cut-off) in the interval 
(A, A + dA) at time t is related to the distribution at the 
initial time to through, nh{A^t) = nh[A-\- X^it — to), to]. 
The two natural choices for the initial states are equilib- 
rium at To = Tc or To ^ oo (all high temperature initial 
conditions become equivalent to the latter). In the for- 
mer case the distribution rih is known analytically and in 
the latter it is very close to the one of critical percolation 
that is also known exactly [30|. Therefore the dynamic 
distributions are known as well [3l|, [32| . More precisely, 
the equilibrium critical Ising hull-enclosed area distribu- 
tion is nh{A^{)) = Ch/A'^ with Ch = l/(87r>/3) and one 
finds 



nh{A,t) 



Ch 



[A^R^{t)Y^ 



(2) 



for t ^ to. Note that this expression exhibits the scal- 
ing form, nh{A^t) = t~'^f{A/t) since the characteristic 
length scales as t^/^. For a quench from infinite temper- 
ature one observes that the distribution takes the critical 
form of random percolation in a few time steps [sil, |32| , 
even though the initial fraction of +1 (or, equivalent ly, 
— 1) spins, 50%, is smaller than the critical density of 
random percolation on a square lattice, pc ^ 0.59. This 
feature was also exploited in (33| to explain the exis- 
tence of percolating blocked asymptotic states in the 
zero-temperature dynamics of this model. 

Unlike in the q = 2 case, nh{A,t) is not analytically 
known for g > 2. To start with, there is no simple relation 
between the area distribution at time t and the one at the 
initial time to and, moreover, we expect the distribution 
to get scrambled in a non-trivial way during the coars- 
ening process. Thus, in a previous work [lH we studied 
this problem with Monte Carlo simulations. First, we 
confirmed that the coarsening process is characterized 
by a growing length that depends on time as [l8|, [l9| 



\g t 



(3) 



by studying different correlation functions and their scal- 
ing properties. Second, we investigated the hull-enclosed 
area distribution. We found that in the cases in which the 
transition is second order (2 < g < 4) and the initial con- 
figuration is critical. To = Tc, the dynamic distribution 
has a power law tail for sufficiently large areas. Quanti- 
tatively, we found that the exponent is 2 independently 
of q (within our numerical accuracy) while the prefac- 
tor depends on q. Interestingly enough, even though the 
individual area rate of change depends on the number 
of sides, the long-range correlations present in the initial 
state are preserved and the system keeps the distribu- 
tion shape during evolution. Assuming that the number 
of sides in the von Neumann-Mullins equation ([1]) can 
be replaced by a constant mean, n{t) (n), and using 
Eq. (j2j) at t = to we proposed 



nh{A,t) 



[A^R^{t)Y^ 



(4) 
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with c)^^ g'-dependent constants. Within our numerical 
accuracy, this form describes weh the huh-enclosed area 
distribution. The results are very different in the case of 
infinite temperature initial conditions. An uncorrelated 
spin configuration representative of equilibrium at Tq 
oo can be mapped onto one of the random percolation 
model. On the square lattice the species density, p = 
is much smaller than the critical random percolation 
limit and the distribution does not become critical nor 
acquire a power-law decay at any time. Instead, it has an 
exponential tail. A similar behavior is observed for g > 4 
even when the initial state is taken from equilibrium at 
Tc. The short-range correlations initially present become 
irrelevant after a finite time and the system loses memory 
of the initial state. See also the results in [s^. 

In the nearest-neighbour Potts model considered here, 
diversely from the so called cellular Potts model in which 
several layers of interacting neighbours are considered, 
the lattice anisotropy is an important ingredient. This 
fact, together with the small number of colors (thus co- 
alescence may be important) deviate the model consid- 
ered here from the ideal grain growth situation. Un- 
der these conditions, it is interesting to see whether the 
areas and perimeter distributions have any resemblance 
with those obtained with different methods or experimen- 
tally. Moreover, the von Neumann- Mullins law, when ap- 
plied to individual grains, also depends on the ideal grain 
growth hypothesis. Thus, in this paper we investigate ad- 
ditional geometric properties of the non-conserved order 
parameter dynamics of the ferromagnetic Potts model 
through Monte Carlo simulations on a square lattice with 
linear size L = 1000. The statistical analysis was per- 
formed using over 1000 samples, enough to reduce the 
errorbars to values that are smaller than the data points. 
Each case was studied by considering an initial equilib- 
rium state at Tq = Tc or an uncorrelated one (Tq oo). 
The critical state (Tq = Tc and 2 < q < 4) was ob- 
tained by running 500 to 4000 Swendsen-Wang algorithm 
steps, while the uncorrelated one was created by ran- 
domly choosing the state of each of the N = L'^ spins 
on the lattice. After reaching equilibrium, the system 
was suddenly quenched to T/ = Tc/2 where we expect 
pinning and nucleation effects not to be present [lH, [26| 
and the subsequent evolution to be a curvature-driven 
process. Time is measured in Monte Carlo units (MCs), 
that is single spin flips such that, statistically, each 
spin is updated once in every step. 

The paper is structured as follows. In Sec.|TTlwe present 
a qualitative discussion of the problem with emphasis 
upon the role played by the initial conditions. In Sec. IIIII 
we discuss the relation between hull-enclosed areas and 
perimeters and their rates of changes. We underline the 
relation between the domain morphology and the char- 
acteristic length R{t) according to the temperature of 
the initial state. In Sect. IIVI we put the von Neumann- 
Mullins equation to the test. In Sec. |V]we show our nu- 
merical results for the time evolution of the hull perime- 
ter distributions. Finally in Sec. IVII we summarize the 



results and we conclude. 



II. INITIAL CONDITIONS 

The initial conditions, either at critical (2 < g < 4) 
or infinite temperature (V g), have very different charac- 
teristics and thus play an important role in the system's 
evolution after the temperature quench, leading to very 
different dynamic structures. When a 2 < g < 4 model 
is in equilibrium at Tc its configuration has long-range 
correlations. In d = 2 the thermodynamic transition also 
corresponds to a percolation transition and one spanning 
cluster is already present in such initial conditions |35l - 
38]. Naturally, such spanning clusters are not counted 
in our analysis since their perimeters would be severely 
under-estimated due to the finite system size. On the 
other hand, when the system presents a first order phase 
transition {q > 4), the correlations are short-ranged at 
To = Tc and the spanning cluster is not present initially 
nor at any time within the timespan (and system size) 
considered in our study. When the initial state is one 
of equilibrium at infinite temperature, Tq ^ oo, the cor- 
relations are absent and the model can be mapped onto 
random percolation. No spanning cluster is present ini- 
tially for g > 2 and none is generated during evolution 
(the case g = 2 is subtly different, see [3ll-[33t). These 
distinctive features are made evident in the snapshots 
shown in Figs. [2] (g = 3) and[3]((7 = 8), for the two initial 
conditions, infinite (left column) and critical temperature 
(right column). 

For (7 = 3, Fig. O the absence of correlations at 
To ^ oo is clear from panel (a) while the long-range 
correlations and a spanning cluster present at Tc can be 
easily visualized in panel (b). In panels (c) and (e) the 
snapshots show the thermal evolution after a quench from 
Tq ^ oo to T/ = Tc/2 at times t = 2^ and 2^° MCs. In 
(d) and (f) the snapshots show the configurations at the 
same times after a quench from Tq = Tc. The perime- 
ters of the domains at t = 2^^ MCs are shown in (g) and 
(h). The right column panels demonstrate that the seeds 
for the large dynamic domains were already present in 
the initial condition and, although these objects do not 
translate during the evolution, they get rid of small, in- 
ternal domains. Indeed, in the explored timespan, basi- 
cally the modifications are that small domains are erased 
and the walls get smoother. A similar conclusion can- 
not be drawn from the left column snapshots as it is 
difficult to identify a seed for the structures seen at the 
latest time from the very disordered initial condition. A 
common feature of both Tq = Tc and Tq ^ oo evolved 
conflgurations is that domains with small areas have an 
approximately round shape. Instead, for large areas the 
pattern is rougher and one finds that these are sometimes 
stretched. A crossover in the morphology of the domains 
is determined by the characteristic length R{t) and it will 
be discussed in Sec. Illli 

Figure [3] shows a sequence of snapshots of the q = S 
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FIG. 2. (Color online.) Snapshots of the q = 3 ferromagnetic 
Potts model on a square lattice with linear size L = 1000 
after a quench, from initial states prepared at To ^ cxo (left 
column) and To = Tc (right column) , to the final temperature 
Tf — Tc/2. Times shown are, from top to bottom, t = 0, 2^ 
and 2^° MCs, respectively. In panels (g) and (h) the interfaces 
in the configurations at the latest time t — 2^° MCs [(e) and 
(f)] are shown, along with some small thermal fluctuations. 



Potts model with both Tq ^ oo (left column) and 
To = Tc (right column) initial conditions. The quench 
is done to Tf = Tc/2, and the configurations at differ- 
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FIG. 3. (Color online.) Sequence of snapshots for the ^ = 8 
Potts model. The left column displays results obtained by 
using an equilibrium initial configuration at To ^ oo (a) and 
its time evolved after the quench [(c) and (e)]. The right 
column shows configurations found by using an equilibrium 
initial state at To = Tc (b) and after the quench [(d) and 
(f)]. From top to bottom the measuring times are t = 0, 2^ 
and 2^° MCs, respectively, and Tf = Tc/2. In panels (g) and 
(h) we show the domains with areas A < (green, or light 
grey), R^{t) < A < 10R^{t) (white) and A > 10 R^{t) (red, 
or dark grey) for the configurations at t = 1024 MCs. 
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ent times, t = 128 MCs (second row) and t = 1024 
MCs (third row) are shown. Since none of the initial 
configurations is critical, thus having either a vanishing 
(To oo) or a finite (Tq = Tc) correlation length, one 
might expect them to be statistically equivalent. How- 
ever, the dynamically evolved configurations are different 
in at least two senses. First, the domains evolved from 
Tc are larger, as can be clearly seen at naked eye. More 
importantly, the shape of the large domains are very dif- 
ferent. Panels (g) and (h) take the configurations in pan- 
els (e) and (f) and color in the same way domains with 
areas A < cR^ (green), cR^ < A < lOi?^ (white) and 
A > lOR^ (red). In this way we highlight the density of 
domains with small, intermediate and large size as com- 
pared to the typical one, i?^, at the measuring time. On 
the one hand, when the initial state is uncorrelated (left 
column panels) few domains with area A > lOR^ sur- 
vive at the latest time (two) and those have a stretched 
form. On the other hand, when the initial state is critical 
(panels in the right column) more domains remain at the 
chosen time (five), they have larger area and are usually 
less stretched. The difference will be made quantitative 
in Sec. IIII Al It is clear from these figures that it is very 
hard to collect good statistics for large areas. 



III. AREAS AND PERIMETERS 

In this section we analyze the relation between hull- 
enclosed areas and hull perimeters. For any sensible def- 
inition of volume and interface, the volume V of a com- 
pact domain with a compact interface S should satisfy 



V 



P 



d-l 



(5) 



with p a linear dimension and d the dimension of space, 
leading to 



(6) 



However, if a domain has a fractal surface or holes (other 
domains) in its interior, Eqs. (j5j) are not necessarily 
valid [5[. In practice, the presence of holes inside do- 
mains depends on the dimension of space and the num- 
ber of ground states of the system. While for d = 2 and 
q = 2 many domains have other internal domains [32| , for 
d = 2 and g > 3 internal domains are much rarer (this 
can be checked from the snapshots in Fig. [3] for q = 8). 
Therefore, domain areas and hull-enclosed areas differ by 
very little for g > 3 [22[. In the rest of this paper we fo- 
cus on hull-enclosed areas and in this Section on their 
relation to the hull perimeters themselves. 

There are many possible definitions of a hull for a sys- 
tem defined on a lattice. We choose to use the one in 
which the hull perimeter is the number of sites that are 
on the external perimeter of a domain, see Fig. (TJ match- 
ing the nomenclature in [5[. 




(a) 



(b) 



(c) 



FIG. 4. (Color online.) Examples of typical domain struc- 
tures obtained in the simulations. In (a) a regular domain 
with area A p^^ and an = 2. The domain in (b) has a 
rougher morphology and a smaller exponent au certifies this 
feature, 1 < < 2. In (c) a domain with very rough and 
stretched morphology, A ^ p] the exponent au reaches it 
minimum value, = 1, in this case. 



A. Scaled scatter plots 

We analyze scatter plots of areas vs. perimeters and 
their scaling properties. The results found in [32] for the 
Ising case are extended in a natural way to the g' = 3, 4 
cases as discussed below. Further subtleties are found in 
models with g > 4 for which none of the initial states 
(To = Tc and Tq ^ oo) are critical. 

As in the Ising case, and independently of the initial 
condition having long or short range correlations, we can 
distinguish two types of domains depending on their re- 
lation to the characteristic length R{t). Hull-enclosed 
structures with area A < R^{t) have a regular form and 
the area-perimeter relation is A ~ as in Eq. (j5j). In- 
stead, domains with area A > R^(t) exhibit a rough sur- 
face and the area-perimeter relation keeps the power law 
form, A ~ p^^ , although with an exponent ah < 2 that 
depends on the initial condition and q. Figure |4] shows 
examples of domains - taken from our simulations - with 
different exponents, a/^ = 2 (a), 1 < Of/i < 2 (b) and 
ah — ^ (c). Once hull-enclosed areas and perimeters 
are measured with respect to the crossover value, itself 
written in terms of the characteristic length R{t)^ we can 
write the scaling relation 



A_ /p 
R^ ^ \R 



where 



= 2 

<2 



for 
for 



A< 
A> 



R\t) 
R\t) 



(7) 



(8) 



Figure [5] shows the rescaled relation ([7]) for g = 3 af- 
ter a quench from Tq oo to Tj = Tc/2 and several 
measuring times (t = 0, 2, . . . , 2^^ MCs). We used the 
characteristic length R{t) ^ {Xq 1)^1'^ obtained from the 
analysis of the equal time correlation C{r^t)^ as in [22j . 
The data collapse confirms that A scales as R^ and p 
as R. We shall see below that this applies to all q and 
to quenches from different initial conditions. While for 
small areas, A < R^(t)^ the exponent can be obtained 
from a fit performed on a conveniently long interval, for 



-| 0"^ - 10^ 10^ 10^ 10^ - 

lO""" 10° 10^ 10^ 10^ 

p/R{t) 

FIG. 5. (Color online.) Scaling relation between hull-enclosed 
areas and hull perimeters at several measuring times for ^ = 3 
after a quench from To ^ oo to T/ = Tc/2. The inset shows 
the raw data at several times from t = to 2^"^ MCs (bottom 
to top) . The extent up to which the law A ^ (red solid line) 
applies increases with time, demonstrating that the domains 
become more regular in the course of evolution. 



A > R^it) the fitting interval is rather short. Neverthe- 
less, the exponent obtained, ah — 1.06, is so far from the 
regular value 2, that we can definitely conclude that a 
change in morphology operates upon the domains at the 
crossover determined by R(t). The domains with small 
areas are regular and A p^. The ones with large ar- 
eas are rough and stretched, the hull perimeters being 
very close to the hull-enclosed area itself, A p. The 
inset shows the raw data for the area A versus the hull 
perimeter p for t = to t = 2^^ MCs (bottom to top). 
This figure confirms that during evolution the domains 
grow and become more regular. The red solid line is a 
guide-to-the-eye and represents A . 

Next, in Fig. [6] we compare the relation between hull- 
enclosed areas and hull perimeters for systems with differ- 
ent values of q quenched from different initial conditions. 
In all cases we find a power law relation, A p^^ ^ with 
varying values of ah- The group of curves labeled (a) 
corresponds to g = 2, 3 and 8 quenched from To = Tc 
to Tf = Tc/2 and the fit yields ah ~ 1.44. Interestingly, 
even the case q = 8 with short-range correlation presents 
an exponent a/i > 1. On the other hand, when the initial 
state is one of infinite temperature, the exponent ah is 
smaller and decreases with increasing number of states. 
In the first case (b) for g = 2 we find ah ^ 1.12. For 
larger values of q it is harder to be conclusive about the 
actual value of the exponent since the number of domains 
with area larger than R^ decreases with increasing q and 
the fitting interval gets smaller. Within our numerical 
accuracy, for g = 3 (c) ah ^ 1.06 and for = 8 (d) 
ah ^ 1.06 but this value is not precise enough to be 
conclusive. 

Summarizing, in all cases the characteristic length R 
grows as t^/^, the structures become more regular, and 



10° 10^ 10^ 10^ 10^ 10^ 

p/R{t) 

FIG. 6. (Color online.) Scaling relation between hull-enclosed 
areas and hull perimeters at t = 2, . . . , 2^° MCs for different 
models (various q values) and initial conditions. The curves 
labeled (a) correspond to ^ = 2, 3, 8 and To = Tc, those la- 
beled (b) are for q — 2 and To ^ cxo, the ones labeled (c) 
have been obtained using ^ = 3 and To ^ oo and, finally, 
the ones carrying the label (d) are for ^ = 8 and To ^ cxo. 
All quenches are to T/ = Tc/2. When the initial state is at 
the critical point, case (a), the relation A ~ p^^ presents a 
similar exponent for all ^, au — 1.44. Instead, after a quench 
from the To ^ oo uncorrelated state the exponent ah de- 
creases from its Q = 2 value a^, — 1-15 to :^ 1.06. Notice 
that for the case q = 8 and To = Tc, case (a), the observed 
exponent is 1.44 because the measure was taken at a rela- 
tively small time, in which the finite correlations present at 
the initial state are still large compared with the characteristic 
lenght. At later times, this exponent will decrease, eventually 
approaching 1.06, as the characteristic lenght grows beyond 
those initial correlations. 

the power law A ^ p'^ for A < R^ is more evident. Con- 
comitantly, less domains have an area larger than R^. 
For To ^ oo the value of the exponent ah characterizing 
the shape of the large hull-enclosed areas decreases with 
increasing consistent with a/,, ^ 1 for g ^ oo. For 
To = Tc the exponent ah takes a higher value though 
still smaller than 2 for large areas and, within our nu- 
merical precision, no dependence on g for g = 2, 3 and 8 
was observed. 



IV. THE VON NEUMANN-MULLINS LAW 

The von Neumann-Mullins law, Eq. ([1]), predicts that 
each hull-enclosed area has a different rate of change de- 
pending only on its number of sides (and not on its area). 
This equation is not fully general; it has been derived 
under certain assumptions that include, as in the Allen- 
Cahn case [27], the fact that the domain wall should be 
close to flat away from the triple points for g > 2, a 
feature that can be achieved at long times and for long 
interfaces only. The law ruling the dynamics of small 
areas with highly curved interfaces or large areas with 
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1000 2000 3000 4000 5000 6000 
A 

FIG. 7. (Color online.) Rate of change of the individual 
hull-enclosed area at time t = 1024 MCs after a quench from 
To ^ oo to Tc/2 {q = 8). The data are grouped by the 
number of sides n = 4 (red, or dark grey) and n = 8 (green, 
or light grey). The vertical dashed line indicates R^{t) oc t. 



very rough walls is not known in general (in the Ising 
case, however, some exact results for small areas are 
known |40| ) . In this Section we numerically analyze 
the validity of the von Neumann-Mullins law by putting 
Eq. ^ and its average over small {Aq < A < cR^) or 
large {A > cR^) areas to the test (Aq is a small area 
cut-off and c a tunable parameter). We estimate the in- 
dividual area time-derivate as ^^(t + 1) — Ai{t). 

In Fig. [71 we examine Eq. ([T]) by plotting dA/dt vs A 
for n = 4 and n = 8 both at t 1024 MCs in a g = 8 
Potts model quenched from To ^ oo to Tc/2 (similar el- 
lipsoidal clouds are obtained for To = Tc and for ^ = 3). 
The vertical dashed line shows R^{t) ex t. The large ver- 
tical spreading of the data points let us conclude that the 
von Neumann-Mullins law does not hold strictly in our 
case, contrary to what was found in a phase field simu- 
lation of ideal grain growth [4l] and in a Potts models 
with very large q and microscopic updates modified to 
ensure a local q = Q constraint and thus accelerate the 
evolution (42| . 

Figure [8] (top panel) shows the area rate of change 
averaged over small areas: 



dA 
~dt 



1 



EdAi 



Ao<Ai<cR'^ 



(9) 



where N^{t) is the number of areas obeying Aq < Ai < 
cR^. When one area leaves the defining interval [Ao, ci?^], 
it is no longer taken into account in the sums; on the 
other hand, new areas can enter this interval. These are 
the reasons for the time-dependence in A^<. Each curve 
in the figure corresponds to different values of q and tem- 
perature of the initial conditions: (a) g = 2 at To ^ oo, 
(h) q = 3 at To = T^, (c) q = 3 at To ^ oo, (d) g = 8 
at To = Tc and (e) (7 = 8 at To ^ oo. In cases (a) 
and (b) the averaged area rate of change has reached a 
constant (within the fluctuations). In the other cases. 
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FIG. 8. (Color online.) Top panel: The area rate of change 
averaged over small areas, Ao = 8 < A < cR^, for (a) ^ = 2 
and To ^ 00 (red), (b) ^ = 3 and To = Tc (blue), (c) ^ = 3 
and To ^ (X) (green), (d) ^ = 8 and To = Tc (pink), and (e) 
q = S and To ^ 00 (cyan). In all cases the averaged rate 
of change seems to asymptotically reach a negative constant. 
Bottom panel: The time-dependent average number of sides, 
(n(t))< for the same parameters is shown with the same color 
and label code. See the text for the discussion. 



although the average continues to evolve within the ex- 
plored time window, the variation is much smaller than 
the total mean area. In all cases {dA/dt) ^ is negative 
and, according to the von Neumann-Mullins's law, this 
implies an average number of sides smaller than 6. The 
fluctuations observed mainly in the case q = 3 are the 
result of coalescence and dissociation processes that be- 
come increasingly rare as the number of ground states of 
the system increases. 

The results in the top panel of Fig. [8] are to be com- 
pared to those in the bottom panel, where the small-area 
average of the time-dependent number of sides, (n(t))<, 
is displayed. The results are qualitatively and quantita- 
tively consistent with the averaged von Neumann-Mullins 
law 



^)^. I ((„(.),, -6) 



(10) 



The values for \q found by comparing the two sets of 
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FIG. 9. (Color online.) The same as in Fig. [8] but averaging 
over large areas, A > ci?^, for the same choice of parameters 
and color code used in that figure. Notice the order of the 
labels in the top panel. 



curves are consistent with Xq decreasing for increasing 
q and independent of the initial temperature To, as is 
also found from the analysis of the space-time correlation 
function [22|. 

The same analysis performed on the large areas is much 
more delicate. To start with, it is difficult to collect good 
statistics since only a few sufficiently large areas remain 
in the samples at long times, especially for large values of 
q. Figure [9] present data for the same set of parameters 
as in Fig. [HI Clearly, the fluctuations in the averaged area 
rate of change are much larger than when the average was 
performed over smah areas. Still, one can argue that the 
rate of change asymptotically approaches a constant in all 
cases. However, and contrary to the small area case, the 
constant is negative only in the cases in which the initial 
state is critical {q = 3 and 4, not shown, at Tq = Tc) and 
for the special case q = 2 with Tq ^ oo. In ah other cases 

= 3, To ^ oo (c), ^ = 8, To = Tc (d) and ^ = 8, To ^ 
oo (e)] the constant is positive. In cases (a) and (b) one 
finds agreement between the negative averaged area rate 
of change and the fact that the averaged number of sides 
of these areas is smaller than 6. Moreover, the Xq values 
found, within our numerical accuracy, are consistent with 
the ones stemming from the analysis of the data averaged 




FIG. 10. (Color online.) Total number of hulls with n sides 
for ^ = 3, To ^ oo (solid symbols, upper group of curves) 
and To = Tc (empty symbols, lower group of curves). The 
dotted lines indicate Nh = 1 (horizontal) and n = 6 (vertical). 
The different symbols correspond to different measuring times 
as given in the key. The distribution has weight well above 
Nh — 1 for n > 6 when the initial conditions have finite 
correlation length. 



over small areas. 

The remaining three cases have to be discussed sepa- 
rately. While the signs are consistent, i.e. positive av- 
eraged area rate of change and average number of sides 
larger than 6, the trend of these curves is not systematic 
and shows that such large structures do not follow the lin- 
ear behavior predicted by the von Neumann-Mullins law. 
This might be due to different reasons. For instance, in 
the case q = 3 with To ^ oo (c) we ascribe this feature to 
the fact that the long-time domain structure is special in 
this case, with many large domains with a rather rough 
surface, as shown in panel (e) of Fig. [2] and confirmed by 
the fact that ah is very close to one, see Fig. O in this 
regime. 

Figure [10] display the total number of hulls with n 
sides, Nhin^t)^ in a log- log scale, showing the behavior 
at four different times, ioi q = 3 after a quench from 
To ^ oo (solid symbols) and To = Tc (empty symbols) 
to Tf = Tc/2. The dotted horizontal and vertical lines 
indicate Nh = 1 and n = 6, respectively. For initial 
conditions with finite correlation length (upper group 
of curves) many interfaces have more than 6 sides and 
the average value (n) is less representative of the fluc- 
tuating behavior. The distribution tail seems exponen- 
tial. Instead, for critical initial conditions (lower group 
of curves), the number of huhs with 6 or more sides be- 
comes negligible during evolution and the distribution 
decays as a power law. Interestingly, for To ^ oo the 
curves are monotonically decreasing for n > 4 while for 
To = Tc they are monotonically decreasing for n > 2 
(note that for the case q = 3 shown here, n only takes 
even values). The former feature is experimentally estab- 
lished although it is not captured by simplified random 
neighbor models [isl, Q • 
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FIG. 11. (Color online.) The full averaged area rate of change 
as a function of the number of sides at t = 64 MCs for the 
parameters given in the key. The solid lines correspond to the 
function f(x) ^ with 7=1, 1.5, 2, from bottom to top, 
and serve as guides to the eyes. 



FIG. 12. (Color online.) Rescaled hull perimeter distribution 
for ^ = 3 at several times (t = 2, 2^, ... , 2^^ MCs), after a 
quench from To = Tc to T/ = Tc/2. The sohd black line 
represents the theoretical prediction, Eq. (|12p . valid for small 
areas, A/R^ < 10 with an = 2 and for large areas, A/R^ > 10, 
with ah ^ 1.44. 



The summary of our analysis of the von Neumann- 
Mullins law is given in Fig. [11] where we averaged the 
area rate of change over all areas with given value of n 
and displayed the results as a function of n for the q = 3 
and q = 8 models quenched from Tq = Tc and Tq ^ 00. 
The data are measured at t = 64 MCs. Although the 
n — 6 dependence is observed for both values of q and 
the different initial conditions, the behavior is linear in 
the averaged data only up to n 10 (lower solid straight 
line). Deviations appear for larger values of n in all cases. 



V. HULL PERIMETER DISTRIBUTIONS 

In a previous work, Ref. [22|, the hull-enclosed area 
distribution, n^(A, t), was studied for the g-state Potts 
model after temperature quenching the system from 
above the critical temperature into the ordered phase. 
We now describe the time evolution of the hull perime- 
ter distribution, nh{p^ t), following this quench and relate 
it to nh{A,t). The form of the distribution depends on 
the initial condition and on the value of q. In particular, 
there is a dependence on the morphological characteris- 
tics of each domain, that are different for structures with 
large or small sizes, when compared to the characteristic 
lenght R{t). 



A. Critical initial condition 

For 2 < g < 4 the transition is continuous and the 
initial state at Tq = Tc is critical; hence all initial dis- 
tributions, i.e. of areas and perimeters, should have 
a power law tail. The functional form of the equilib- 
rium number density of hull-enclosed areas at t = is 
nh{A,0) ((7-l)4^V^^ El- Since one and only one hull 



is associated to each hull-enclosed area, the number den- 
sity of hull perimeters should be linked to the one of hull- 
enclosed areas according to nh{p^ 0)dp = nh{A^ 0)dA. Us- 
ing A ~ p"^, a relation that is valid one-to-one, one finds 



(11) 



where 7^ = ah{q — and C/i = <^/i + 1- The exponent 

is the equivalent to the Fisher exponent r for the area 
distribution now describing the perimeter distribution. 
In two dimensions is linked to the fractal dimension 
of the perimeter as [5] C/i = 1 + '^/Dh- In the special 
case g = 2 at Tc one has Dh = 11/8 [35]. 

Using the fact that areas and perimeters are related 
one to one also dynamically, nh{p-,t)dp = nh{A^t)dA^ 
and the power-law relation between areas and perime- 
ters, Eq. ([7]), the time-dependent number density of hull 
perimeters can be easily written as 



ill) 



ah-l 



(12) 



The value of ah depends on which of the two regimes, 
p < R and p > R^ one focuses on, with a/^ = 2 in the 
former and a/^ < 2 in the latter. A similar argument 
was used in [32] and [45| to obtain the perimeter number 
density in the 2d Ising model with non-conserved and 
conserved order parameter dynamics, respectively. 

Figure [12] shows the hull perimeter number density, 
scaled by the characteristic length R{t)^ for the critical 
model with q = 3 after a quench to Tf = Tc/2. The long- 
range correlations present in the critical initial state are 
preserved during evolution and lead to the power law tail 
in the distribution. The exponent of this tail is (h ^ 2.5 
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FIG. 13. (Color online.) Hull perimeter distribution for q — 8 
at several times given in the key after a quench from equilib- 
rium at To ^ oo to T/ = Tc/2. The declivity of the envelope, 
Ch 3.0, is a direct consequence of dynamical scaling. 



in agreement with the theoretical estimate C/i = + 
1 c:^ 1.44+1, see Fig. [6l For small areas, p <C i^, one 
finds E?nh{p-,t) — ^qp/R also in good agreement with 
the analytic prediction. Thus, equation ([12]) describes 
well these two limits. 



B. Non-critical initial condition 

Being far from critical percolation, the initial state at 
To ^ oo for a system with q> 2 does not become critical 
during its evolution and the area (and consequently the 
perimeter) distribution does not develop a power law tail. 
In Fig. [13] we use log-log scale to show the q = 8 hull 
perimeter distribution after a quench from Tq ^ oo to 
Tf = Tc/2 at several times. Similarly to what happens 
with the hull enclosed area distribution for g > 2 [22| . 
there is ap~^ envelope that is a consequence of dynamical 
scaling. 

For very large after the quench, domains start to 
increase from localized density fluctuations that are ran- 
domly scattered throughout the system^Jn a similar way 
to the Avrami-Johnson-Mehl method |4g] used to pro- 
duce Voronoi diagrams. Thus, neglecting coalescence ef- 
fects that are only important for small it can be con- 
jectured that the set of domain centers (defined, for ex- 
ample, as the center of mass of a domain) may be ap- 
proximately described as a set of random points from 
which Voronoi domains are drawed. Indeed, as we show 
in Fig. [H for g = 8 (see also Ref. [22[), the distribu- 
tion of areas after a temperature quench from Tq ^ oo 
to Tc/2, is well described by the generalized Gamma (or 
Voronoi) distribution [46|, widely used in grain growth 
literature EHMl: 



FIG. 14. (Color online.) Hull-enclosed area distribution for 
q = 8 at several times after a quench from equilibrium at 
To ^ cxo to T/ = Tc/2, rescaled by R^(t). The solid line is a 
fit with Eq. ^ and a 1.21, b ^ 1.22, c 0.6. Inset: the 
same data in a linear scale, showing that for small areas, the 
collapse improves at large times. 



with a 1.21, b :^ 1.22, c 0.6. Simpler versions of 
the above function have also been proposed, with two 
(c = 1) or a single (c = 1 and a = b) parameter. Indeed, 
the values that we get for a and b are almost the same, 
but c 7^ 1. 

Once the distribution of areas is known, in princi- 
ple one could relate it, as in the preceeding section, to 
the perimeter distribution using Eq. ([7]), obtaining, once 
again, a Gamma distribution: 



n.(.,t)-i?(|)"^""exp[-6(|) 



(14) 



exp 



(13) 



Differently from the area distribution, the dependence on 
ah characterizes two distinct regimes, = 2 for p < R 
and ah < 2 for p > R. As one can see from Fig. [TS] the 
data are well described by the above distribution. How- 
ever, the fit parameters are not only different from those 
obtained from the area distribution but differ as well be- 
tween the regions with small and large p. One possible 
reason is that although Eq. ([T3]) is a good approximation, 
the dynamics of a growing domain is not equivalent to the 
Avrami-Johnson-Mehl method. If the Voronoi cells of the 
domain centers are drawn, there will be a large superpo- 
sition with the actual domains, but those regions close to 
the interfaces could be wrongly assigned because, differ- 
ently from Voronoi diagrams in which the borders are as 
straight as they can be on a lattice, a coarsened domain 
has rougher interfaces. Thus, the deviations are on the 
scale of the perimeter lenght and, although they may not 
be noticeable for area distributions, they become visible 
in the perimeter distributions. 
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FIG. 15. (Color online.) Hull perimeter distribution for q — 8 
at several times after a quench from equilibrium at To ^ oo 
to Tf = Tc/2, rescaled by R{t). The solid line is Eq. (O 
with ah = 2 ior p < R and ah — 1.06 for p > R. The 
fit parameters, however, besides being different from those of 
Fig-im differ also for p < R and p > R. Inset: the same data 
in a linear scale. 



VI. CONCLUSIONS 

We studied several geometric properties of the nearest- 
neighbour g-state Potts model quenched below its phase 
transition. In particular we analyzed the von Neumann- 
Mullins law, the hull perimeter length distribution and 
the relation between hull perimeters and their enclosed 
areas. Diversely from the so called cellular Potts model, 
in which several layers of interacting neighbours are con- 
sidered in order to decrease the lattice anisotropy, here it 
is an important ingredient. Moreover, coalescence is also 
relevant when q is not too large, as the values considered 
here. These two factors, necessary for the single domain 
von Neumann-Mullins law, are here violated and deviate 
the model from the ideal grain growth conditions. 

After having confirmed that all dynamic observables 
satisfy dynamic scaling with respect to the growing 
length R{t) ^ (A^t)^/^, we studied the fractal proper- 
ties of the hulls by plotting A/R'^{t) against {p/R{t))^^, 
Eq. [71 We found that for p/R{t) <C 1 interfaces are 
smooth with ah = 2 while for p/R{t) 1 they are fractal 
with ah being smaller for To ^ oo than for To = Tc and 
decreasing with q. Thus, the behavior of objects with lin- 
ear sizes that are smaller that the typical growing length 



have a rather different dynamics, morphology and statis- 
tical properties than those whose linear sizes are larger 
than the dynamical length. 

The linear proportionality between dA/dt and n (in- 
dependently of A), the von Neumann-Mullins law, does 
not hold for each individual area for the nearest neigh- 
bour Potts model in the temperature interval considered 
here after the quench. This may be due to thermal fluctu- 
ations occuring along the perimeter that may mask the 
curvature driven contribution. We then examined the 
law on the mean, by averaging over areas that are smaller 
or larger than the typical one, R^{t). In the small area 
regime the linear von Neumann-Mullins law is verified on 
the mean while for larger areas it is not, with deviations 
from linearity depending on the value of q and initial con- 
ditions. The same separation of regimes is found when 
averaging over all areas that have the same number of 
sides. The behavior at large values of n seems to be cap- 
tured by a modified power law involving ah but our data 
are not extensive enough to test this conjecture. 

Finally, we studied the hull perimeter distribution dis- 
tinguishing critical (To = Tc for q> ^) from non-critical 
(all other cases, especially large q for To oo). In the 
former case we found that the scaling function can be de- 
rived with an argument that combines the use of Eq. ([7]) 
with the form of the distribution of hull-enclosed areas 
found in Ref. [22|. For non-critical initial conditions we 



found, instead, that the Gamma function commonly used 
in studies of ideal grain growth describes well both the 
hull-enclosed area and perimeter distributions, although 
the expected relation between the exponents of both dis- 
tributions is not obeyed. The reason is due to the fact 
that the perimeter distribution is much more sensitive 
to the differences between the actual domains and their 
Voronoi approximations, as discussed in section IVBI 

There are, however, several statistical properties of in- 
terface sides and topological properties for cellular sys- 
tems that were not yet fully studied for the nearest- 
neighbour Potts model, and will be the subject of a future 
publication. 
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